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ABSTRACT 

Context. The stable element hafnium (Hf) and the radioactive element thorium (Th) were recently suggested as a suitable pair for 
radioactive dating of stars. The applicability of this elemental pair needs to be established for stellar spectroscopy. 
Aims. We aim at a spectroscopic determination of the abundance of Hf and Th in the solar photosphere based on a CO s BOLD 3D 
hydrodynamical model atmosphere. We put this into a wider context by investigating 3D abundance corrections for a set of G- and 
F-type dwarfs. 

Methods. High-resolution, high signal-to-noise solar spectra were compared to line synthesis calculations performed on a solar 
CO s BOLD model. For the other atmospheres, we compared synthetic spectra of C0 5 BOLD 3D and associated ID models. 
Results. For Hf we find a photospheric abundance A(Hf)=0.87 ± 0.04, in good agreement with a previous analysis, based on ID 
model atmospheres. The weak Thn 401.9 nm line constitutes the only Th abundance indicator available in the solar spectrum. It 
lies in the red wing of an Ni-Fe blend exhibiting a non-negligible convective asymmetry. Accounting for the asymmetry-related 
additional absorption, we obtain A(Th)=0.09 ± 0.03, consistent with the meteoritic abundance, and about 0.1 dex lower than ob- 
tained in previous photospheric abundance determinations. 

Conclusions. Only for the second time, to our knowledge, has a non-negligible effect of convective line asymmetries on an abun- 
dance derivation been highlighted. Three-dimensional hydrodynamical simulations should be employed to measure Th abundances 
in dwarfs if similar blending is present, as in the solar case. In contrast, 3D effects on Hf abundances are small in G- to mid F-type 
dwarfs and sub-giants, and ID model atmospheres can be conveniently used. 

Key words. Sun: abundances - Stars: abundances - Hydrodynamics 



1. Introduction 

The determination of the ages of the older stellar populations in 
the Galaxy is fundamental for constraining its buildup mecha- 
nisms and determining the associated timescales. A variety of 
dating techniques have been developed, the majority relying 
on comparing the colour-magnitude-diagram of a stellar pop- 
ulation with theoretical isochrones or "fiducial lines" derived 
from other observed populations (typically globular clusters). 

Long lived unstable isotopes of heavy elements such as 
238 U and 232 Th provide an independent method of deter- 
mining stellar ages. Such elements are generally produced 
through neutron capture (specifically rapid neutron capture or 
r-process) and then decay. Their abundance at the present time 
thus provides a direct estimate of the time passed since they 
were produced, if their initial abundance can somehow be in- 
ferred. This is typically accomplished by measuring some other 
stable element with a common nucleosynthetic origin together 
with the radioactive element. Europium (Eu) has long been a 



typical choice, since its abundance relative to the abundance of 
thori um can be established from r-prqcess enrichment model s 
(seelFrancois et al1ll993l ICowan et alJI 19991 llvans et aT] |2006 ) . 
ISneden et alJ (I2000h used the [Th/Eu] abundance ratio to es- 
timate the age of the old, metal-poor globular cluster M15 to 
14 ±3Gyr. The study was made easier by the fact that Ml 5 
shows strong r-process element e nhancement ([Eu/Fe]=1.15, 
rTh/Fel=0.81. ISneden et"aDl2000h . Nevertheless, Eu presents 
some drawbacks as a cosmochronometric reference element: 
while easily measurable, its atomic mass is substantially differ- 
ent from the one of Th and U. This generally increases uncer- 
tainties. In fact, decaying and stable reference elements should 
both be formed in the same event, which - for heavy n-capture 
species - can be safely assumed only if the atomic masses 
of the two elements are similar. This means in practise that, 
when a stable reference element is used, it should be as mas- 
sive as possible. Alternatively, where U is measured, the U/Th 
ratio can be used since the two have very different half-life 
times. The measurement of U is nevertheless extremely diffi- 
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Fig. 1. Four different observed solar spectra, displayed for four 
different spectral windows, each containing a Hfn line, indi- 
cated by a thin dashed vertical line. For clarity, the intensity 
spectra, labelled as DI (Delbouille Intensity) and NI (Neckel 
Intensity), are shifted down by -0.2 units, while the disk- 
integrated spectra, KF (Kurucz Flux) and NF (Neckel Flux), 
are shifted up by 0.2 units for the 356.1 nm and 391.8 nm lines. 
Except for the DI spectrum, the original normalisation has been 
retained. Note that the different spectra do not always agree 
perfectly. 
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Fig. 2. The Fe i+Ni i blend in the 3D simulation (solid line) is 
compared to the IDlhd profile (dash-dotted line), scaled in or- 
der to have the same central residual intensity as the 3D profile. 
The difference of the two spectra is shown as well as a Th 3D 
line profile. 
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Fig. 3. The 3D profile of (Fei+Nii) + (Th n+Co i+ V n) blend 
in solid black is over-imposed on the IDlhd fit (green/grey 
dashed line). In the plot is also over-imposed (black dashed 
line) the IDlhd with the same abundance as the 3D profile as 
well as the 3D Th n+Co i+V n profile (dotted line). 



cult a nd has been possible only in two metal-poor r-enhan ced 
stars (ICayrel et al.l(200lblHill et aIll2002l:lFrebel et al.ll2007h . 

Recently, hafnium (Hf, Z=72) has been suggested as a suit- 
able stable reference element. The Hf production rate via the 
r-process is believed to be strictly tied to the one of thorium, so 
that their abundance ratio is almost constant regardless of the 
neutron density in the generating Super Novae (SN) (Kratz et 
al. l2007h . Hafnium is more massive than Eu (stable Eu isotopes 
are 15 1 and 153, while Anf ~ 178) and can be measured from a 
set of lines in the blue part of the visual spectrum (the most suit- 
able ones being at 391.8 and 409.3 nm) for which accurate tran - 
sition data have been recently measured (lLawler et a l. 2007). 
On the down-side, Hf is believed to be efficiently produced by 
s-process (about 55% of the meteoritic Hf content in the solar 
system is believed to have been produced by the s-process, see 
lArlandini et al.|[l999h . As a consequence, it can be effectively 
employed as a reference only where s-process enrichment can 
be considered negligible, or if the s-process contribution can be 
reliably modelled. 

Although in principle promising, Hf/Th dating is not free 
of drawbacks. In the first place an un certainty of 0.1 dex in 
[Th/Hf] implies about 4.7 Gyr (see also lLudwig et al. Il2008ah 
uncertainty in the dating. Second, the theoretical prediction of 
strong ties between Th and Hf production needs to be verified 
by observations. Third, the applicability of the method to nor- 
mal (i. e. not r-process enhanced) stars needs some testing: es- 
pecially Th lines become relatively weak and significant blends 
can be prejudicial especially in more metal-rich objects. 

With the prospect of using hafnium and thorium for dating 
purposes, we decided to determine the solar hafnium and tho- 
rium abundances using hydrodynamical model atmospheres. 
We anticipate that 3D effects are found to be non-negligible 
for Th, and mainly due to the convection-induced asymmetry 
of blending lines. 

In the solar spectrum, the measurable lines of Hf are from 
the singly ionised atom, and are often blended with other 
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Table 1. Hf n lines considered in this work. 



Table 2. The CO' BOLD models considered in this work. 
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photospheric absorption lines. As a consequence, solar pho- 
to spheric_j&undance determinations are scarce in the litera- 
ture. Russell ( 1929) determined t he photospheric haf nium solar 
abundance to be 0.90. Recently. lLawler et alj d2007l) measured 
Hf lifetimes for several lines, four of which were used to deter- 
mine a solar ph o tosphe ric hafnium abundance of 0.88 + 0.08. 
lAndersen et al.l (1 19761) con sidered nine lines o f Hf n (two of 
which are in common with lLawler et al.1 12007) and obtained 
from a Kitt Peak National Observatory centre-disc spectrum a 
solar photospheric hafnium abundance of 0.88 + 0.08. Values 
of the meteori tic hafnium abundance ci ted in literature include: 
0.73 ± O. Ol.lAnders & Grevessd(ll989l): 75 + 0.02 Grevesse 
& Sauval ( 1998); and 0.77 + 0.04|L^ddersJ d2003|) — 

Thorium is measurable in the Sun from one single line 
of Thu at 401.9 nm, which is heavily b lended. So l ar tho rium 
abundances in literature are v ery rare. iHolwegerl (1 19801) de- 
rived a Th abundance of 0.16. lAnders & Grevessd (1 19891) re- 
port abundances of 0. 12 +0.06 from the 401.9 nm line d ue to an 
unpublished measurement of Grevesse. lLoddersI d2003l) recom- 
mends the meteoritic value A(Th)=0.09 + 0.04 from Grevesse 
et al. (119961). The meteoritic value is 0.08 + 0.02 ( Anders & 
Grevesse l 19891) . inferred from the mass-spectroscopicallv mea- 
sured r atio of 0.0329 Th atoms to 10 6 Si atoms ( Rocholl & 
Jochum 1993) in carbonaceous chondrites. 



2. Atomic data 

W e considered the sam e four Hf n lines with the same log gf 
as lLawler et al.l ((2007) have been considered (see Table [TJ. 
Line parameters f or the 401.9 nm Thn and blending lines in 
this range are from ldel Peloso et"aTI d2005al) . 



3. Models 

Our analysis is based on a 3D model atmosp here computed 
w ith the CQ 5 BOLD code dFrevtag et al.ll2002t Wedemever et 
al. 120041) . In addition to the CQ 5 BOLD hydrodynamical sim- 
ulation we use several ID mo dels. More de tails about the ap- 
plied mod els can be found i n Caffau et al. (20 071) and Caffau 
& Ludwig H2007I) . As in previous works we call (3D) the ID 
model derived by a horizontal and temporal averaging of the 
3D C0 5 BOLD model. The reference ID model for the com- 
putation of the 3D abundance corrections is a hydrostatic ID 
model atmosphere computed with the LHD code (see Caffau 
& Ludwig 120071) . hereafter lDi »n model. lDi »n employs 
the same micro-physics (opacities, equation-of-state, radiative 
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Note: Cols.(l), (2), and (3) state the atmospheric parameters of the 
models; col(4) the number of snapshots N t considered; col(5) the 
time interval covered by the selected snapshots; col(6) a characteristic 
timescale t c of the evolution of the granular flow. 



transfer scheme) as C0 5 BOLD, and treats the convective en- 
ergy transport with mixing-length theory. Here, we set the 
mixing-length parameter to 1.5. We also use the solar ATLAS 9 
model comput ed by Fiorella Caste ll{_] and the Holweger-M uller 
solar model dHolwegerlll967l;lHolweger & Mullerlll974l) . The 
spectrum synthesis code we employ is Linfor3E0 for all mod- 
els. 

The 3D C0 5 BOLD solar model covers a time interval of 
4300 s, represented by 19 snapshots. This time interval should 
be compared to a characteristic time scale related to convec- 
tion. Here we choose the sound crossing time over one pressure 
scale height at the surface f c = Hp/c (Hp: pressure scale height 
at T ross = 1, c: sound speed) which amounts to 17.8 s in the so- 
lar case. Note, that this time scale is not the time over which the 
convective pattern changes its appearance significantly. This 
time sc ale is significantly longer than f c ; e.g JWedemever et alj 
(2004) obtained from simulations a autocorrelation life time of 
about 120 s for the convection pattern. However, while f c does 
not measure the life time of the convective pattern as such, it 
provides a reasonably accur ate scaling of its life time with at- 
mospheric parameters (e.g JSvensson & Ludwig 2005). Hence, 
it can be used for the inter-comparison of different 3D models. 

The particular solar model used in this paper has a higher 
wavelength resolution (12 opacity bins) with respect to the 
model we used t o study the solar phosphorus abundance 
(CaffaueFaD|2007). The basic characteristics of the 3D models 
we consider in this work can be foun d in Table [2] 

As in our study of sulphur (see Caff au & Ludwig 120071) 
and phosphorus (see lCaffau et al.ll2007t) . we define as "3D cor- 
rection" the difference in the abundance derived from the 3D 
C0 5 BOLD model and the IDlhd model, both synthesised with 
Linfor3D. We consider also the difference in the abundance de- 



1 http://wwwuser.oats.inaf.it/castelli/sun/ap00t5777g44377klasp.dat 

2 http://www.aip.de/ mst/Linfor3D/linfor_3D_manual.pdf 
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rived from the 3D C0 5 BOLD model and from the (3D) model. 
Since the 3D C0 5 BOLD and (3D) models have, by construc- 
tion, the same mean temperature structure, this allows us to 
single out the effects due to the horizontal temperature fluctua- 
tions. 

4. Observational data 

The observational data we use consists of two high-resolution, 
high signal -to-noise ratio spectra of centre disc solar inten- 
sity: that of iNeckel & La bs ( 1984) (hereafter referred to as the 
"Neckel intensity spectrum") and that of (Delbouille, Roland 
& Neven Tl973h ( hereafter referred to as the "Delbouille inten- 
sity spectrum"E| ). We further use t he disc-integrated solar flux 
sp ectrum of Necke l & Labsl ( I 19841) and the solar flux spectrum 
of iKuraczl (12005 fl referred to as the "Kurucz flux spectrum". 

5. Data analysis 

5.1. Hafnium 

The computation of a 3D synthetic spectrum with Linfor3D is 
rather time consuming and the spectral-synthesis code can han- 
dle easily only a few tens of lines. The log gf values of the lines 
blending the Hf lines and lying nearby are not well known, 
and the synthetic spectra do not reproduce the observed solar 
spectrum well in that region. To fit the line profiles Lawler et 
al. (I2007I) changed the log gf of some lines blending the Hf n 
lines, but did not make available the values used. In any case, 
solar log gf values derived from ID and 3D fitting are not guar- 
anteed to be the same. Our computational power is insufficient 
at the moment to enable us to fit 3D synthetic spectra to an ob- 
served spectrum by changing the log gf values of various lines. 
Consequently, we could not produce solar log g f for the blend- 
ing lines and fit the line profiles with a 3D model grid. Instead 
of fitting astrophysical solar log gf of the blending features, 
we introduced blending components in the process of measur- 
ing EW. We thus took advantage of the deblending mode of the 
IRAF task splot. The corresponding Hf abundance is then de- 
rived from the curves of grow th computed with L infor3D for 
the various lines. The results of Lawl er et al.l (I2007I) are used as 
a check for our measurements. The abundances we obtain from 
the Holweger-Miiller (HM) model are in agreement with their 
measurement. The results for all spectra and individual lines 
are presented in Table [3] 

We determine the following hafnium abundances: 0.870 + 
0.038 if we consider only the disc-centre intensity spectra, 
0.818 + 0.005 according to the solar flux spectra, and 0.856 + 
0.040 if we consider both the intensity and the flux spectra. As 
final uncertainty we took the line to line RMS scatter over all 
spectra. Equivalent widths can be measured in both flux spec- 
tra only for the 409.3 nm hafnium line, which is nearly free of 
blends. Only the Neckel flux spectrum can be considered for 
the Hf line at 39 1 . 8 nm, due to the fact that the larger broaden- 
ing of flux spectra with respect to intensities makes the mea- 
surement more difficult, if not impossible. Additionally, NLTE 

3 http://bass2000.obspm.fr/solar_spect.php 

4 http://kurucz.harvard.edu/sun.html 



corrections may be different in intensity and flux spectra. For 
these re asons, and to have a consistent comparison with the re- 
sults of lLawler et al.l d2007), we determine the solar hafnium 
abundance from the average value stemming from the disc- 
centre s pectra. Our re sults are in good agreement with the re- 
sults of (Lawler et ail d2007l) : we obtain A(Hf)=0.878 + 0.040 
from the Holweger-Miiller model, considering only disc-centre 
spectra, which is to compare to their value of 0.88 + 0.02. 

The 3D abundance corrections for Hf are very small. The 
highest value, 0.037 dex, is the one for the strongest line. As 
previously mentioned, a 3D correction is the difference in 
abundance obtained from the 3D model and the IDlhd model. 
If one chooses an ATLAS or the Holweger-Miiller model as a 
ID reference model, the 3D correction would be even smaller 
or negative, respectively. 

The (3D) and IDlhd solar model temperature structures 
are similar. The relative contribution to the 3D correction due 
to different average temperature profiles is small. The 3D cor- 
rection is larger for stronger lines, quite insensitive to micro- 
turbulence, of the same order of magnitude for flux and inten- 
sity, and slightly higher for intensity. 

5.2. Thorium 

Singly ionised thorium is the dominant ionisation stage in the 
solar atmosphere, and thorium is measurab le in the Sun only 
through the 401.9 nm Thn resonance line. Holweger (1980) 
detected also the 408.6 nm Thn line, however in our view this 
line is too weak and blended to be used for quantitative analy- 
sis. The 401.9 nm line happens to fall on the wing of a strong 
blend of neutral iron and nickel. Additionally, it is blended with 
a Co i line (about 25% of the EW of the Th n line) and a weaker 
V ii line (1/5 of the Co line, 1/9 of the Th line in EW). The Co i 
and V ii lines are very close in wavelength to the Th n line, and 
it is therefore practically impossible to disentangle the blend. 
Consequently, we do not measure the contributi on in EW due 
to tho rium, but instead use the EWs measured bv lLawler et alj 
(I 19901) to obtain A(Th) (se e Tabled first two lines). However, 
while lLawler et all d 19901) considered only the contribution of 
Co i to the blend, we consider also the contribution of V n. To 
do so we computed from the CO s BOLD solar m odel the EW of 
the Co i +V ii blend, using the atomic data from del Peloso et al. 
(I2005al) . and subtracted this value from the EW of lLawler et all 
(Il990l) . In this manner, we obtained a smaller solar A(Th) (see 
last two lines of Table|4]i. The 3D corrections listed in Table|4] 
are deduced from the EWs alone, and we shall call them in the 
following "intrinsic" 3D corrections. These are the corrections 
that would be found if the Th n line were isolated. Since it is 
instead found in the middle of a rather complex spectral region, 
other 3D-related effects arise, which are discussed later. 

The last two columns of Table [4] provide estimates of the 
statistical uncertainty of the theoretical EW predicted by the 
3D model and the associated uncertainty in the hafnium abun- 
dance. The 3D model provides only a statistical realisation of 
the s tellar atmosphere, a nd is thus subject to limited statistics 
(see lLudwig et al. 120 08b). However, as necessary for a reliable 
derivation of abundances these uncertainties are insignificant 
and only added for completeness. 
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Fig. 4. Fits of the synthetic spectra to the observed solar spectra. On the left hand side, the synthetic spectra have been computed 
from the CO s BOLD model, on the right hand side from the corresponding IDlhd model. On either side, from top to bottom, the 
spectra are: the Delbouille disc centre intensity, the Neckel disc centre intensity, the Neckel flux and the Kurucz flux. On each 
panel the sub-components used to synthesise the feature are shown individually: the light dotted line is the Th n+Co i+V n blend, 
the darker dotted line is the Ni i+Fe i blend, the black dashed lines are, from blue to red, the Mn i+Fe i blend the artificial Fe n 
line and the Co 1 401 .93 nm lines. 
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Table 3. Solar hafnium abundances from the various observed spectra. 



spec wavelength EW A(Hf) <x Ref. 3D-1D LH d °"<ew> °~a 

nm pm 3D (3D) IDlhd IDatlas IDhm dex dex % dex 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



Nl 


338.9 


0.489 


0.91 1 


0.895 


0.883 


0.900 


0.914 


0.02 




0.029 


0.29 


0.0013 


DI 


338.9 


0.477 


0.899 


0.883 


0.871 


0.886 


0.903 


0.02 


0.91 


0.029 






Nl 


356.2 


0.845 


0.819 


0.801 


0.781 


0.799 


0.821 


0.01 




0.037 


0.31 


0.0014 


Dl 


356.2 


0.851 


0.822 


0.805 


0.785 


0.803 


0.824 


0.02 


0.85 


0.037 






NF 


391.8 


0.318 


0.823 


0.827 


0.815 


0.836 


0.848 


0.05 




0.008 


0.26 


0.0011 


Nl 


391.8 


0.302 


0.911 


0.899 


0.893 


0.912 


0.925 


0.02 




0.018 


0.27 


0.0012 


Dl 


391.8 


0.274 


0.867 


0.854 


0.849 


0.868 


0.881 


0.05 


0.91 


0.018 






KF 


409.3 


0.317 


0.814 


0.818 


0.805 


0.826 


0.838 


0.07 




0.010 


0.26 


0.0011 


NF 


409.3 


0.318 


0.816 


0.820 


0.806 


0.828 


0.840 


0.07 




0.010 






Nl 


409.3 


0.292 


0.886 


0.874 


0.865 


0.884 


0.899 


0.02 




0.020 


0.26 


0.0011 


Dl 


409.3 


0.268 


0.847 


0.835 


0.827 


0.846 


0.860 


0.04 


0.86 


0.020 







Note: Abundances and their uncertainties - if not noted otherwise - are given on the usual logarithmic spectroscopic scale where A(H)=12. 
Col. (1) is an identification flag, DI means Delbouille and Nl Neckel disc-centre spectrum, NF Neckel and KF Kurucz flux spectrum. Col. (4): 
hafnium abundance, A(Hf), according to the CO 5 BOLD 3D mod el. Cols. (5)-(8): A(H f) from ID models, withf, micro-turbulence, of 1.0 km s~ ] 
. Col. (9): uncertainties related to EW measurement. Col. (lO) lLawler et al.l d2007t) . Col. (11): 3D corrections. Col. (12): uncertainty of the 
theoretical equivalent width due to the limited statistics obtained from the 3D model. Col. (13): uncertainty in abundance due to the model 
uncertainty in equivalent width. 

Table 4. Solar thorium abundances derived from the EWs of lLawler et al.l (Il990h . 



spec 


EW 




A(Th) 




a 


3D-(3D) 


3D- 1D LHD 


C <EW> 


OA 




pm 


3D 


(3D) 


IDlhd 


1D H m 








% 


dex 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Int 


0.415 


0.168 


0.164 


0.148 


0.190 


0.022 


0.004 


0.0207 


0.33 


0.0015 


Flux 


0.505 


0.135 


0.151 


0.130 


0.172 




-0.015 


0.0052 


0.34 


0.0015 


Int* 


0.355 


0.095 


0.092 


0.075 


0.118 


0.026 


0.003 


0.0199 






Flux* 


0.430 


0.060 


0.074 


0.054 


0.096 




-0.015 


0.0060 







Notes: The contribution to the EW of Co and V, according to the 3D simulation, has been been from the EWs in the last two rows (labelled 
with *). (1) identification flag to distinguish disc-centre to flux spectrum results; (3) thorium abundance, A(Th), according to the CO 5 BOLD 
3D model; (4)-(6) A(Th) from ID models, with^, micro-turbulence, of 1.0km s -1 ; (7) the uncertainty due to the EW; (8)-(9) 3D corrections; 
(10) statistical uncertainty of the theoretical EW; (11) corresponding uncertainty in abundance. 

Table 5. Solar thorium abundance derived from line profile fitting. 







3D 






LHD 








A(Th) 


SHIFT 


v br 


A(Th) 


SHIFT 


v br 


3D-1D LHD 




dex 


km s _I 


km s -1 


dex 


km s~' 


km s _I 


dex 


Int. (Neckel) 


0.088 


-0.05 


1.49 


0.160 


-0.53 


2.77 


-0.720 


Int. (Delbouille) 


0.080 


+0.01 


1.86 


0.168 


-0.47 


2.99 


-0.088 


Flux (Neckel) 


0.056 


-0.39 


1.68 


0.168 


-0.47 


4.11 


-0.112 


Flux (Kurucz) 


0.067 


-0.38 


2.82 


0.175 


-0.46 


4.71 


-0.108 



Note: (2)-(4) results of the fitting: A(Th), shift, and broadening using a 3D grid; (5)-(7) same for IDlhd grid; (8) gives the difference in 
abundance A(Th) 3D -A(Th) 1DLHD . 



Since the Thu line lies on the red wing of the stronger 
Fe i+Ni i blend, one can expect that the line asymmetry of the 
latter contributes in a non-negligible manner to the measured 
EW. Up to now it has generally been considered that convective 
line asymmetries have no impact on the derived abundances, 



which is true, if the line is isolated. ICavrel et alj d2007l) have 
demonstrated that such asymmetries are non-negligible for the 
derivation of the 6 Li/ 7 Li isotopic ratio. From the purely mor- 
phological point of view, the case of the Thu resonance line 
is similar to that of the 6 Li i resonance line: a weak line which 
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falls on the red wing of a stronger line, which is asymmetric 
due to the effects of convection. In order to investigate this ef- 
fect one has to consider the full line profile and resort to spec- 
trum synthesis. 

To gain some insight into the problem we began by compar- 
ing synthetic profiles, which included only the Fe i+Ni i blend 
and the Th n+Co i+V n blend, computed from the C0 5 BOLD 
and IDlhd models. In Fig|2] the IDlhd flux profile of the 
Fe i+Ni i blend has been scaled (since the line strength is dif- 
ferent in ID and 3D) in order to have the same central depth 
of the 3D profile. The asymmetry of the red wing is obvious. 
The difference between the 3D and the IDlhd profiles is shown 
in the dotted line of the figure. The 3D thorium profile is also 
shown as a dashed line and it makes clear that the line asym- 
metry does indeed contribute to the measured EW. The EW 
of this difference is 0.134 pm, which translates into a differ- 
ence in abundance of +0.13 dex for flux. For intensity, A(EW) 
= 0.088pm implies AA(Th)=+0.10dex in Th abundance. 

On a level of slightly higher sophistication with respect to 
the above naive argument, we fitted the 3D synthetic profile 
with a grid of IDlhd synthetic spectra of differing Th abun- 
dance. We scaled the IDlhd Fe+Ni blend in order to induce 
it to have the same central residual intensity as the 3D pro- 
file. One of the results of the fitting is shown in Fig. [3] the 3D 
profile (black solid line) with A(Th)=0.09 when fitted with a 
IDlhd grid (dashed green/grey line) gives a higher abundance 
A(Th)=0.171. The use of a (3D) grid leads to an even higher 
thorium abundance: A(Th)=0.227. In the plot the IDlhd flux 
profile, with the same abundance as the 3D synthetic spectrum, 
is shown, as well as the 3D thorium profile. For the disc-centre 
intensity, with the same procedure as used for flux, the fitting 
of a 3D profile assuming a thorium abundance of 0.09 with 
a IDlhd grid implies A(Th)=0.110, while with a (3D) grid 
A(Th)=0.155. The difference of the result of the fit when us- 
ing (3D) or IDlhd grids is due to the superposition of the ef- 
fect of the line asymmetry and the "intrinsic" 3D corrections 
( A(Th) 3D- A(Th) m) . 

The Fei+Nii line asymmetry and the 3D corrections act 
in opposite directions: the line asymmetry produces a higher 
A(Th) from ID analysis, while the 3D correction is positive, 
thus A(Th) from a ID analysis is smaller. However, the 3D 
correction in intensity is between 3 and 4 times larger than for 
flux (see Table|4|i. 

In Fig. [4] we show the results of fitting the whole spectral 
region both with 3D (left hand side) and IDlhd (right hand 
side) profiles. The synthetic profiles have been built by adding 
5 components: a Mnl+Fel blend centred at 401.89 nm, the 
Fe i+Ni i blend, the Th n+Co i+ V n blend, an artific i al Fe n line 
at 401.9206, initial ly suggested by iFrancqis et al] (|T993) and 
adopted also by Idel Peloso et al.1 (l2005ah . and the two Co i 
lines at 401.93. Of these components the Thu+Co i+V n blend 
was computed for three different Th abundances, all the other 
components were instead computed for a single abundance and 
were scaled by the fitting procedure in order to best reproduce 
the observed spectrum. 

The results of the fitting are summarised in Table [5] From 
the table one sees that the total 3D correction, taking into ac- 
count both line asymmetry and "intrinsic" 3D correction is 
from -0.07 to -0.09 dex for the intensity spectra and -0. 1 1 dex 



for the flux spectra. This difference is expected, given that the 
"intrinsic" 3D correction, as given in Table HI is about +0.02 
dex for intensity and about +0.005 dex for flux. 

It is quite interesting to note that in the case of Th the line 
asymmetry implies a total 3D effect which is not only much 
larger, but also in a different direction, from what would have 
been implied by the "intrinsic" 3D correction. To our knowl- 
edge this is the second time such an effect has been highlighted. 
These "total" 3D corrections are only slightly larger than what 
was estimated above from the fit of the 3D Fe i+Ni i blend plus 
Thn+Coi+Vn blend, with IDlhd profiles. This suggests that 
the use of the simpler procedure is still capable of providing 
the correct order of magnitude for the "total" 3D correction. 
Finally we point out that, as can be seen from the left hand 
side of FiglU it is obvious that the IDlhd profiles are unable 
to provide the correct shape for the Th feature. The shape is re- 
produced in a much more satisfactory way by the 3D profiles. 

We conclude this section by investigating the iron line at 
401.9nm, which is the stronger component of the Fei+Nii 
blend. The 3D synthetic line shows a strong asymmetry; the 
difference 3D- ID is very pronounced on the red wing if we 
make an appropriate shift and broadening to the ID profile so 
as to over-impose the blue wings of the 3D and ID synthetic 
spectra. While the flux of the 3D synthetic spectrum is shifted 
to the blue by O.lkms -1 , the intensity spectra, correspond- 
ing to different inclination angles, show a shift ranging from 
-0.5 km s" 1 to +0.15 km s~' . 



6. The 3D effects for other models. 

6.1. Hafnium 

Beyond the solar model as such we investigated the hydrody- 
namical effects on some other 3D models of dwarf and sub- 
giant stars (see Table |2|. We restricted our investigation to 
solar or -1.0 metallicity, because for lower metallicity these 
hafnium lines are no longer measurable in dwarfs. Even the 
detection at -1.0 metallicity is only possible for the strongest 
of the hafnium lines, but the two metallicities can be used for 
interpolation to intermediate metallicity. We find that 3D and 
ID abundances are generally very similar, making the 3D cor- 
rections insignificant for all models and all lines we consider. 

Hf ii is the dominant species of hafnium in the photosphere 
of these models. Both the 5900 K models predict that at least 
98% of the Hf is ionised. The hotter models foresee an even 
lower fraction of Hf i. The coolest models at 5000 K predict 
that at least 75% of the Hf is ionised. The results are similar 
when considering the (3D) or the IDlhd temperature profile 
for the atmosphere. 

If we compare the horizontal average temperature profile of 
the 3D model to the IDlhd model we can see that they are very 
similar. There is no pronounced cooling of the 3D atmospheres 
relative to IDlhd in radiative equilibrium for the same atmo- 
spheric parameters which is known to be prominent in more 
metal-poor models. 
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6.2. Thorium 

As for hafnium and as in the solar case, the "intrinsic" 
3D corrections are small for the other models. The hydro- 
dynamical effects are instead due to the line asymmetry, 
like in the Sun. When we fit, as we do for the Sun, a 
3D flux profile, with meteoritic or scaled meteoritic A(Th) 
for the metal-poor models, with a IDlhd grid we obtain 
for A(Th) a value higher by 0.08-0. 12 dex considering the 
models: 5500K/3.50/0.0, 5500K/3.50/-1.0 5900K/4.00/0.0, 
5900K/4.00/-1.0, 5900K/4.50/0.0, and 5900K/4.50/-1.0 (see 
Tabled . One has to keep in mind that such differences translate 
into systematic age differences of almost 5 Gyr in the context 
of nucleocosmochronology. 

7. Conclusions 

Our determination of the photospheric hafnium abundance, 
A(Hf)=0.87 + 0.04, derived from a 3D hydrodynamical 
CO s BOLD solar model, i s very close to the abundance deter- 
mination of ILawler et al.1 (120071) obtained using the Holweger- 
Miiller solar model. The EWs of the hafnium lines are not 
very sensitive to the hydrodynamical effects. A difference of 
the same amount of the 3D correction can be found if we 
compare the A(Hf) obtained using two different ID models, 
such as ATLAS and Holweger-Miiller model or IDlhd and 
the Holweger-Miiller model. These difference are related to 
different opacities and physics used to compute the models. 
Synthetic profiles of hafnium lines computed from 3D models 
show a pronounced asymmetry, but this effect is not relevant 
for abundance determination as long as only the EW is consid- 
ered. 

For thorium, as for hafnium, 3D effects on the EWs are less 
than or comparable to the uncertainties related to the EW mea- 
surement, so that they are negligible when determining A(Th). 
However, for the solar photospheric thorium determination the 
asymmetry of the strong Fei+Nii line blending the 401.9 nm 
Th ii weak line implies a non-negligible effect on the derived 
Th abundance and is in fact the dominant source of the "total" 
3D effect. The asymmetry is a hydrodynamical effect, and we 
find a difference in the A(Th) determination of about -0. 1 dex 
when considering this asymmetry. Our results for the Sun and 
for other 3D dwarfs and sub-giants models imply that Th can- 
not be reliably measured without making use of hydrodynami- 
cal simulation if the Th line is weak and blended. A reappraisal 
of existing Th measurements in the light of hydrodynamical 
simulations is warranted for dwarfs and sub-giants. 
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Wavelength 


EW 




A(Hf) 




3D-<3D> 


3D-1Dlhd 


f <EW> 


0~A 


nm 


pm 


3D 


<3D> 


1D LHD 






% 


dex 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


5060/4.5/0.0 


















339.0 


0.718 


0.870 


0.842 


0.829 


0.028 


0.041 


0.05 


0.0002 


356.2 


1.550 


0.872 


0.809 


0.796 


0.062 


0.076 


0.09 


0.0003 


391.8 


0.405 


0.870 


0.854 


0.841 


0.016 


0.029 


0.07 


0.0003 


409.3 


0.410 


0.870 


0.855 


0.841 


0.016 


0.029 


0.07 


0.0003 


5060/4.5/- 1.0 


















339.0 


0.151 


0.870 


0.866 


0.852 


0.004 


0.018 


0.11 


0.0005 


356.2 


0.395 


0.870 


0.864 


0.852 


0.007 


0.018 


0.12 


0.0005 


391.8 


0.080 


0.870 


0.870 


0.855 


0.001 


0.015 


0.11 


0.0005 


409.3 


0.081 


0.868 


0.867 


0.853 


0.001 


0.015 


0.11 


0.0005 


5870/4.5/0.0 


















339.0 


0.551 


0.870 


0.862 


0.858 


0.008 


0.011 


0.20 


0.0008 


356.2 


1.150 


0.870 


0.857 


0.849 


0.012 


0.020 


0.22 


0.0010 


391.8 


0.330 


0.871 


0.867 


0.869 


0.004 


0.002 


0.18 


0.0008 


409.3 


0.336 


0.870 


0.867 


0.867 


0.003 


0.003 


0.18 


0.0008 


5920/4.5/-1.0 


















3390. 


0.095 


0.871 


0.872 


0.868 


-0.001 


0.003 


0.57 


0.0025 


3562. 


0.230 


0.870 


0.880 


0.881 


-0.010 


-0.011 


0.67 


0.0029 


3918. 


0.053 


0.869 


0.872 


0.873 


-0.003 


-0.004 


0.55 


0.0024 


4093. 


0.054 


0.866 


0.869 


0.870 


-0.003 


-0.004 


0.56 


0.0026 


6230/4.5/0.0 


















339.0 


0.444 


0.870 


0.865 


0.859 


0.005 


0.012 


0.22 


0.0009 


356.2 


0.916 


0.870 


0.865 


0.855 


0.004 


0.015 


0.25 


0.0011 


391.8 


0.275 


0.870 


0.866 


0.870 


0.004 


-0.001 


0.19 


0.0008 


409.3 


0.281 


0.870 


0.866 


0.869 


0.003 


0.001 


0.20 


0.0009 


6240/4.5/-1.0 


















339.0 


0.073 


0.869 


0.876 


0.870 


-0.007 


-0.001 


0.19 


0.0008 


356.2 


0.173 


0.871 


0.889 


0.889 


-0.018 


-0.018 


0.20 


0.0009 


391.8 


0.042 


0.866 


0.874 


0.876 


-0.008 


-0.010 


0.17 


0.0007 


409.3 


0.044 


0.875 


0.884 


0.885 


-0.009 


-0.010 


0.17 


0.0007 


6460/4.5/0.0 


















339.0 


0.374 


0.871 


0.869 


0.864 


0.002 


0.006 


0.23 


0.0010 


356.2 


0.764 


0.870 


0.871 


0.863 


0.000 


0.007 


0.26 


0.0011 


391.8 


0.239 


0.869 


0.868 


0.877 


0.002 


-0.008 


0.20 


0.0009 


409.3 


0.245 


0.871 


0.869 


0.877 


0.002 


-0.006 


0.20 


0.0009 


6460/4.5/-1.0 


















339.0 


0.059 


0.867 


0.883 


0.875 


-0.016 


-0.007 


0.25 


0.0011 


356.2 


0.137 


0.868 


0.897 


0.897 


-0.029 


-0.029 


0.27 


0.0012 


391.8 


0.035 


0.865 


0.881 


0.882 


-0.016 


-0.017 


0.21 


0.0009 


409.3 


0.036 


0.866 


0.883 


0.883 


-0.016 


-0.017 


0.21 


0.0009 


5500/3.5/0.0 


















339.0 


1.500 


0.868 


0.855 


0.850 


0.014 


0.019 


0.23 


0.0010 


356.2 


2.780 


0.871 


0.845 


0.830 


0.026 


0.041 


0.29 


0.0012 


391.8 


0.987 


0.870 


0.865 


0.867 


0.005 


0.003 


0.21 


0.0009 


409.3 


1.000 


0.868 


0.864 


0.863 


0.004 


0.005 


0.22 


0.0009 


5500/3.5/-1.0 


















339.0 


3.160 


0.870 


0.870 


0.863 


0.001 


0.007 


0.29 


0.0013 


356.2 


0.758 


0.870 


0.874 


0.862 


-0.003 


0.008 


0.34 


0.0015 


391.8 


0.182 


0.871 


0.874 


0.872 


-0.003 


-0.001 


0.28 


0.0012 


409.3 


0.186 


0.871 


0.875 


0.870 


-0.004 


0.001 


0.29 


0.0013 



Note: The atmospheric parameters are given as r eff /logg/[M/H]. cols. (3)-(5) are the A(Hf), according to CO s BOLD 3D, (3D) and 1D LH d 
model; cols. (6)-(7) are 3D corrections, 3D-(3D) and 3D-1D L hd respectively; col. (8) statistical uncertainty of the theoretical EW; col. (9) 
corresponding uncertainty in abundance. 



